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ABSTRACT  -  The  Lee-Nikodem  equations  of  motion  of  elastic  plates  are  solved 
for  the  case  of  vibrations  of  an  AT-cut  quartz  strip,  with  free  faces  and 
edges,  at  frequencies  up  to  and  including  the  third  harmonic  thickness-shear 
overtone. 

\ 

1.  Introduction 

About  30  years  ago,  A.W.  Warner  [1]  developed  a  high  precision 
crystal-plate  resonator  utilizing  ‘he  third  harmonic  overtone  of  thickness- 
shear  vibration,  i.e.  a  mode  involving  a  thickness-shear  motion  with  three 
nodes  across  the  thickness  of  the  plate  rather  than  the  one  node  of  the 
fundamental  thickness-shear  mode.  At  about  the  same  time,  equations  were 
developed  which  extended  the  classical  (Lagrange-Germa in-Cauchy)  range  of 
frequencies  to  include  that  of  the  fundamental  thickness-shear  mode ;  but  it 
was  not  until  much  later  that  Lee  and  Nikodem  [2,3]  formulated  equations 
suitable  for  studying  vibrations  at  frequencies  of  the  harmonic  overtone 
modes  of  thickness-shear. 

In  the  present  paper,  the  Lee-Nikodem  third-order  equations  are 
solved  for  a  case  of  rotated-Y-cut  quartz  plates  with  free  faces  and  a  pair 
of  parallel,  free  edges.  The  results  of  computations  for  the  AT-cut  plate 
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are  presented  for  vibrations  in  the  neighborhood  of  the  frequency  of  the 
fundamental  thickness-shear  mode  and  in  the  neighborhood  of  the  third 
harmonic  overtone.  The  differences  between  the  two  exhibit  some  of  the 
reasons  for  the  higher  stability  of  the  latter. 


2.  Lee-Nikodem  Equations 


To  obtain  two-dimensional  equations  of  motion  of  plates  from  the 
three-dimensional  equations  of  linear  elasticity,  Lee  and  Nikodem  start  with 
an  expansion  of  the  three-dimensional,  rectangular  components  of  displacement, 
u-  ,  j*l,2,3,  in  series  of  trigonometric  functions  of  the  thickness-coordi- 

J 

nate,  x2  »  of  the  plate: 


uH  *  u<">  cos  n6  , 

J  n*0  J 


A  C  C  8  S  5  i  C  ■ 

I'll  ™  ( ; "■ 

mo  jt 

'■  i".''  :.;icu  . 


1A 


(1) 


where  the  ujn^  are  Independent  of  X£  and 


6  «  n(l  -  x2/b)/2 


L. 


fl 


(2) 


i _ 


in  which  b  is  the  half-thickness  of  the  plate.  The  functions  cos  n$  give 
the  shapes  of  the  simple  thickness-modes  of  an  infinite.  Isotropic  plate  with 
free  faces  at  x2  *  +  b  . 

The  expression  (1),  for  the  u.  ,  Is  substituted  in  the  variational 

J 

equation  of  motion  [4]: 
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•4(T  1j,1  •  o»j>«“jdv  *  0 


(3) 


where  the  are  the  components  of  stress,  p  is  the  mass  density  and 
V  is  the  volume.  The  integration  is  performed  over  the  thickness  of  the 
plate  and  leads  to  stress-equations  of  motion  of  order  n  ;  which  are, 
omitting  the  terms  accounting  for  surface  tractions. 


-  (mr/2b)  =  enpU<n)  ,  (4) 

where 

Ti")  '  !>-'  /  T,j  cos  ne  dx2  ,  T<J)  •  b-'  /  Ttj  sin  ne  d*., . 
-b  -b 

(5) 

and  en  *  2  for  n  *0  and  e„  *  1  for  n  >  0  .  (Corrections  of  [3]  by  a 
factor  of  2,  for  n*0,  were  kindly  supplied  by  Professor  Lee). 

The  three-dimensional  strain-displacement  relations, 

Sij  *  (uj,i  +  u1,d)/2  * 


become,  with  (1), 


S 


cos  ng  + 


sin  n6)  , 


(7) 


where 
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3.  Reduction  to  a  Special  Case 


The  example  to  be  studied  is  one  of  steady  vibrations  at  frequencies 
high  enough  to  include  the  third  harmonic  overtone  of  the  thickness-shear 
family  of  modes  of  an  AT-cut  quartz  plate  bounded  by  free  faces  at  x2  =  +  b 
and  free  edges  at  x^  +  a.  The  modes  are  to  be  straight-crested  along  x3 
and  antisymmetric  with  respect  to  both  Xj  and  x2  .  Thus,  we  take,  of  (1), 
only 

u-|  *  (u^  cos  6  +  u<3>  cos  36)  e^w*  , 

u2  “  (u2°^  +  u22^  C0S  e1Wt  *  O4) 

u3  *  (u(0)  +  tt(2)  cos  26)  eiwt  , 

in  which  the  ujn^  depend  only  on  x^  .  The  second  term  in  u^  accommodates 
the  third  harmonic  overtone  thickness-shear  mode. 

What  remain  of  the  stress -equations  of  motion  (4)  are 


+  2 


P*2u<°> 


0  , 


2pA<°> 


0  , 


T12J  “  ^/b)T22)  +  p0,2u22)  =  0  * 
Tll!l  “  (37r/2b)T21)+  P“2ul3^  *  0  • 


(15) 


and  the  only  non-zero  components  of  strain  are,  from  (8)  and  (14), 
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c(°)  = 

u(°) 

U3,l  » 

c(°)  - 

^6 

u(°) 

u2,l  * 

s(l)  . 

bl 

uO) 

ul,l  * 

56”  * 

(ir/2b)u^1) 

II 

CM 

— '  m 

Ul 

u(2) 
u3,l  • 

* 

u(2) 
u2,l  * 

■ 

(ir/b)u^2)  , 

542>  ■ 

(n/b)u^  , 

c(3)  . 

u(3) 
ul,l  » 

5<3)  . 

(3rr/2b)uj3* 

Nine  constants  of  elasticity,  referred  to  axes  in  and  normal  to  the 
plane  of  the  plate  (with  an  axis  of  two- fold  symmetry  of  the  elastic  pro¬ 
perties  of  quartz),  enter  into  the  present  example.  As  computed  by  A.  Ballato 

2  -9 

[5]  from  R.  Bechmann's  [6]  principal  constants,  they  are  (in  N/m  *10  ): 


11  ’ 

86.74 

C12  * 

-8.260543013 

U 

in 

in 

o 

68.80698505 

H 

CM 

CM 

129.7663387 

II 

CM 

U 

5.700423178 

c66  * 

29.01301496 

44  * 

38.61152627 

O 

N 

-3.654869573 

C56  * 

2.533571817 

Of  the  remaining  twelve  constants,  four  (Cj3>  c23,  c33,  c43)  do  not  enter 
into  the  present  example,  as  the  modes  are  Independent  of  x3  ;  and  the  others 

^c15*  C25*  C35*  C45*  C16*  C26*  C36*  C46^  are  zer0  *or  the  rotated-Y"cuts  oY 
quartz. 

Lee  and  Nlkodem  Introduce  a  low  frequency  correction  factor  k-j  and  a 
high  frequency  correction  factor  k2  .  The  former  appears  as  a  factor  of  A^q 
in  the  strain  energy  density  appropriate  to  the  present  example: 
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2U  = 


2^55S50,S50)  *  ‘66S<0,S<°>  *  ^6S(0)S<»),  ♦  c„S<’>S<’> 
+  C55S52)S52)  +  c66S62’S62)  +  ^AA'A  +  cnS{3)s{3> 


bb  b  b  bb  b  b  bb  b  ~b  111  1 

+  c  5<2>5<2>+c  5(2)S(2)*2c  5<2>5<2>  *  c  stDjO)  +  c  5(3)g(3) 
c22s2  *2  c44s4  4  ‘c24a2  s4  66  6  H  c66s6  6 

+  2klA10<c66S60)  +  c56S50)>56,)  *  2A12<C66S62>  +  CS«SS2),5«,) 

+  2A2)(c,25<2)  4  c14S(2))S,(,)  *  2A23(c,25<2)  4  c,45<2))S<3) 

+  ^CUSA  +  c56S50)>563)  *  2A32<c66S62)  *  C56S52))563)  • 

(1 


where 


10  *  4/ir  * 

A12  *  "4^3ir  * 

Ag-j  *  8/3ir 

>23  “  “8/5ir  , 

A30  *  4^3ir  *  . 

A32  =  12/5ir 

08) 


The  correction  factor  kg  ,  In  the  present  example.  Is  Inserted  as  a  divisor 
of  the  term  2pw2U2^  In  the  first  of  (15). 

Adjusted  values  of  k-j  and  k2  ,  as  supplied  by  Professor  Lee,  are 


*2/8  » 


*  0.901, 


(19) 


From  (12),  (17)  and  (16),  the  surviving  stress-displacement  relations 


are 
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fO) 

T12 

T(2) 

J22 

T0) 

'23 

t(3) 

12 


ZhA°\  +  cs6(u<“]  ♦  *  b-’u<3) 

^i-66(u<»j+kl.-'u(')t  6-'u») 

cllul!l  +  (8/3bHc-|2U22^  +  c14u32^  ’ 
C55u3,l  +  c56^u2,l  "  (2/3b)uj^  +  (18/5b)u 
C56U3J  +  C66^U2J  ”  (2/3b)u{^  +  (W5b)u 


]  , 


]. 


(3) 

1 

(3) 

1 


3. 

3. 


(20) 


cnu{  j  +  (8/5b)(c^2u2  +  c14u32^  » 

( 4ki /tt) ( c56u^®^  +CggU2j)  +  (ir/2b)Cg6uj^  -  (4/3Tr)(c5gU^2]  +c66u2j)  » 

(8/7r)c12(u|Jj/3  -  uj3j/5)  +  (ir/b)(c22u^2^  +  c24u^)  , 

(8/tt)cu(u1(1j)/3  -  uj  j/5)  +  (Tr/b)(c24i42)  +  c44u^2))  , 

(4/3n)(c5gU^|  +  c66u<°>)  +  (12/5ir)(CggU^2|  +  CggU^2])  +  (3Tr/2b)Cgguj3^ 


The  displacement  equations  of  motion,  to  be  solved,  are  obtained  by 
substituting  the  stress-displacement  relations  (20)  Into  the  stress-equations 
of  motion  (15)— with  k2  Inserted  in  the  first  of  (15),  as  mentioned 
previously. 

Finally,  the  edge  conditions  are 


•(0)  _  T(0)  „  T(l) 


'13 


12 


'll 


t(2)  .  t(2)  .  t(3)  . 
13  "  r12  "  'll 


*  0 


on 


x-|  =  +  a 


(21) 


^  ,  «  i 
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4.  Dispersion  Relation 


In  (14)  we  take,  omitting  the  factor  e1w  , 


,(°) 

2 

=  A^0)  sin  SX,  , 

u<2> 

-  a!2^ 
«2 

sin  £xj 

'l0) 

*  A^  sin  £x^  , 

42) 

=  a(2) 

A3 

sin  Cx1 

'S11 

■  a|^  cos  Sx-j  , 

“!3) 

s  A0) 
A1 

cos  Cx1 

and  substitute  first  In  (20)  and  the  result  in  (15)  to  produce  a  set  of  six 
simultaneous,  homogeneous,  linear  algebraic  equations  in  the  six  amplitudes 
A<">  of  (22): 

J 

a11A2°^  +  a12A3°^  +  a13Al^  +  0  +  0  +  a16Al'3)  *  0 

a12A2°)  +  a22A3°)  +  *zA']  +  0  +  0  +  a26Al3)  *  0 

a13A2°^  +  a23A3°^  +  a33Al  +  a34A2^  +  a35A32)  +  0  *  0  ^ , 

0  +  0  +  a34A1^  +  a44A2^  +  a45A3^  +  a46Al^  =  0 

0  +  0  +  a35A|^  +  a45A^  +  a55A^2)  +  ag6Aj3^  «  0 

a16A20)  +  a26A30)  +  0  *  a46A22)  +  a56A32)  +  a66A{3)  *  0 


The  coefficients  ,  made  dimensionless  and  real  by  some  manipulations 
of  the  equations,  are 
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3'j'j  =  2(z  -a  /kg) *  *  2CggZ  ,  3^  *  4k-|Z  /tt  ,  a^g  -  4z  /it  , 

a22  *  2{c55z2-fi2),  323  *  4k^CggZ  /tt,  3gg  =  AiEggZ  /it, 

a33  *  (c11z2  +  l  -a2)z2,  334  *  -4(1+4c12)z2/37T»  835  *  4(4c14  -  c56)z2/3tt  , 

(24) 

®44  -  z2  +  4c22-ft2,  345  =  4cg4  +  c56z2,  34g  *  4(9  +  4c12)z2/5tt  , 

a55  =  E55z2  +  ^44  -  fl2  •  a56  =  4(4E14  +  9£56)z2/5Tr* 

a66  =  (Ellz2  +  9  '  fl2)z2  • 

where 

✓ 

z  *  2£b/ir  ,  a  =  w/w  *  u2  =  TT2Cg6/4pb2  ,  cpq  =  c^/Cgg  ; 

i.e.  z  is  the  ratio  of  the  thickness  of  the  plate  to  the  half-wave-length 
along  the  plate  and  a  is  the  ratio  of  the  frequency  to  that  of  the  funda¬ 
mental  thickness-shear  mode  of  the  infinite  plate. 

The  determinant  of  the  coefficients  of  the  ,  set  equal  to  zero: 

l*Ml  *  0.  (25) 

in  which  ai4  *  a15  •  a24  *  a25  *  a36  «  0  ,  apq  *  aqp  ,  produces  the  dis¬ 
persion  relation  a  vs.  z  :  a  sextic,  algebraic  equation  in  z2  .  The 
equation  is  the  same  as  (43)  of  [3]  except  for  the  factors  2  in  a^ ,  a12,a22 
as  already  noted  above  in  connection  with  (4).  Also,  here,  all  the  elements 

a  are  real  as  a  result  of  multiplication  of  the  third  and  sixth  rows  and 
P<1 

columns  by  z  . 
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The  six  branches  of  the  dispersion  relation,  computed  on  the  HP-85 
micro-computer,  are  illustrated  in  Fig.  1.  The  characters  of  the  branches 
are  indicated  by  their  identifying  symbols: 

F  =  Flexure 

FS  =  Face-shear 

1^  *  1st  Thickness-shear  (in  x^ -direction) 

2^  =  2ncl  Thickness-shear  (in  x^-direction) 

3-j  =  3  (Harmonic)  Thickness-shear  (in  x^ -direction) 

2.£  ~  2nt*  Thickness-stretch  (in  Xg-di  recti  on) 

The  subscripts  in  the  symbols  1^,  2 3.  3j,  2£  designate  the  direction  of 
displacement  (or  predominant  displacement)  at  z  =  0,  whereas  the  numbers 
themselves  give  the  number  of  nodes  between  Xg  =  +  b.  Thus:  in  2  3  the 
displacement  at  z  =  0  is  predominantly  in  the  direction  of  x^  with  two 
nodes  across  the  thickness  of  the  plate.  Note  that  the  roots  z  for  branches 
F  and  FS  are  real  for  all  Q  ,  but  the  roots  for  the  remaining  four 
branches  may  be  real  or  imaginary,  depending  on  the  frequency.  If  imaginary, 
the  variation  of  displacement  along  is  exponential  or  hyperbolic  rather 
than  trigonometric. 

The  zigzags  in  the  curves  in  Fig.  1  result  from  the  spacing  of  dots 

on  the  cathode  ray  tube  display  of  the  HP-85.  The  figure  is  the  HP-85's  hard 

copy  of  the  CRT  display.  The  roots  z  were  actually  computed  to  an  accuracy 
-9 

of  10  —  a  precision  required  for  their  subsequent  use  in  solving  (34)  and 

(37).  Intervals  of  0.02  in  0  were  employed  for  Fig.  1,  resulting  in  a  com¬ 
putation  time,  for  the  range  0<fl<4  ,  of  about  6  hours  or  about  18  seconds 
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per  root.  The  secant  iterative  method  was  used,  with  starting  values  given 
by  the  following  approximate  formulas,  followed  by  increments  of  10"®  in 


F  :  zf) 

1  .  z2( 

'l  *  z3' 

=  6.42258(1+G)a2[l +(1+K)^]/tt2  , 

(26) 

G 

•  ir2/12(Eu -c22/c22),  K  =  4G(a"2-l)/(l+G)2 

(27) 

FS  :  z2 

=  0.44119  S22 

(28) 

2  •  z2 

z3  ’  z4 

j2.229(a2/a2  -  l),  a  <  a4  , 

'  0.42395(R2/a4  -  1 ) ,  a  >  fi4  , 

(29) 

2  •  z2 

z2  •  z6 

*  I6(a2/flg  -  l),  a  >  ag , 

(30) 

2  2 
°4’«6 

*  2{c22  +  C^  +  [(c22  -  c44)2  +  4c24]^} 

(31) 

3  •  z2 
J1  *  z5 

|0.33799(a2  -  9),  a  <3, 

1 0.40651  (a2  -  9),  a  >3. 

(32) 

These  trial  roots  match  closely  or  exactly  the  roots  of  the  sextic  at  z  =  0 
and  at  a* 0,3  (except  Zg  at  ft* 3)  resulting  in  trial  values  adequate 
for  convergence  of  the  Iteration  for  all  0<a<4  . 
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5.  Frequency  Spectrun 
2 

For  each  of  the  roots  zn  of  (25),  five  amplitude  ratios,  say 


A(0)yA(l).a  A(0)/A(1)-a  A(2>/A(1)«a  A(2)/A(1)»a  A(3)/A(1)-a 

A2  /A1  aln*  A3  /a1  a2n*  A2.  /A1  a3n*  A3  /A1  a4n*  A1  /A1  a5n* 


(33) 


may  be  found  from  five  of  the  six  equations  (23).  Thus,  with  the  third  of  (23) 
omitted,  we  may  write 


all 

a12 


(Vln>  * 

a12(zn°2n* 

+  0 

+  0 

+  a16a5n  * 

<Vln>  + 

a22^zna2n^ 

+  0 

+  0 

+  a26a5n  * 

0  + 

0 

*  *44(zn“3n) 

*  *«S<lna4n> 

+  ,46“5r.  ■ 

0  + 

0 

+  *«5<V‘3n> 

+  *55*lna4n' 

*  *56a5n  * 

~a13 

’a23 

”a34 

“a35 


(34) 


lB^n^ln7  '  "26VV 


This  form  Is  chosen  because  the  zn<*ln,  *n<»2n»  zna3n»  zna4n  and  a5n  are 

real  for  all  fi  ,  as  are  also  the  a  --  as  arranged  previously. 

PM 

With  the  six  zn  from  (25)  and  the  thirty  opn  determined  from  (34), 
we  may  now  write,  in  place  of  (22): 
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.(0) 


•5 


An“ln  s1n  Vi  * 


b 

ui(,)  *  ZJ  a„  cos  Vi  ■ 

n- 1 


.(2) 


SA"‘ 


a4n  s1n  Vi  • 


4°’  '  ZfVzn  s1n  Vi  • 


42’  '  li  V3n  s<"  Vl  •  <35> 

n- 1 


,(3) 


SA"“* 


cos  5„x, 


Upon  substituting  the  displacements  (35)  in  the  formulas  (20)  for  the 
stresses  and  the  results  in  the  edge  conditions  (21),  we  have  the  six  equations 


fiv 


0  ,  m=l . ..6  , 


(36) 


where 


'in 


B2n 


(£562ralo  +  £55zna2n  +  *  2Vl5l/,,)cOS  V 

<Vl„  +  £S62n“2n  *  V  +  2%/w,cos  V  ■ 


*-1 


*>3n  ■  (_£ll2n  +  16ln°3r./3  *  ,6£14zna4n/3)zn  s,n  V  ’ 

b4n  *  *‘4W3”  +  c56zna3n  +  c55Zna4n  *  36c56a5n|,5x^c0s  V  * 

b5n  *  <*4/3"  +  zna3n  +  £S6zn“4n  +  3fa5n/5*)cos  5„l  . 

b6n  ’  (,b£122nV5"  +  16£14zn“4n/5"  +  £Hz2<15r,lz^  s1n  V  * 


in  which  zR  ■  irzn/2  ■  £nb,  A-a/b  and  the  b^  are  real  for  all  Q  . 
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The  roots  l  of  the  equation  obtained  by  setting  the  6x6  determinant 
of  the  coefficients  of  the  AR  ,  in  (36),  equal  to  zero: 


(37) 


produce  the  data  for  plotting  a  frequency  spectrum  fl  vs.  a/b  . 
The  results  of  computations  in  the  two  ranges: 


and 


0.99  <  fl  <  1.01  ,  16  <  a/b  <  24 
2.995  <  fl  <  3.005,  18  <  a/b  <  22 


are  illustrated  in  Figs.  2  and  3.  To  construct  these  figures,  the  six  roots 
zn  of  the  sextic  (25)  were  first  computed  for  a  given  fl.  Then  the  five 
linear  equations  (34)  were  solved  for  the  <Xpn  for  each  of  the  six  zn  and 
the  resulting  combinations  of  Opn  and  zn  substituted  in  the  transcendental 
equation  (37),  after  which  the  range  of  l[-  a/b)  was  traversed  in  steps  of 
0.1  for  Fig.  2  and  0.025  for  Fig.  3  and  the  value  of  Ib^J  computed  at 
each  step.  A  change  of  sign  of  Ib^J  Indicated  a  straddled  root  l  which 

_3 

was  then  determined  to  10  by  successive  linear  Interpolations.  The  process 

•  -5 

was  then  repeated  at  Intervals  of  fl  of  5x10  .  Figs.  2  and  3  required 

about  58  and  49  hours  of  computation,  respectively,  on  the  HP-85. 

In  Fig.  2: 

F22...30  are  overtones  of  flexure 

FS11...15  are  overtones  of  face-shear 

1-|1  Is  the  1st  thickness-shear  (fundamental). 
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In  Fig.  3: 


F62...74 

are  overtones  of  flexure 

FS37...41 

are  overtones  of  face-shear 

1-|33...35 

are  anharmonic  overtones  of  the 
thickness-shear 

1st 

(fundamental) 

2323 _ 27 

are  anharmonic  overtones  of  the 
thickness-shear 

2nd 

transverse 

3^1  is  the 

3™  harmonic  thickness-shear  overtone. 

The  nunbers  following  the  symbols  F,  FS,  1j,  23  and.  31  designate  both  the 
order  of  the  overtone  and  the  approximate  number  of  half-wave- lengths  between 

X1  =  ±  a  • 

Fig.  2  illustrates  the  well  known  phenomenon  of  strong  coupling  of  the 
Is*  thickness-shear  fundamental  with  flexure  overtones  and  weak  coupling  with 
face-shear  overtones.  Fig.  3  shows  that  the  3*^  harmonic  thickness-shear 
mode  has  moderately  strong  coupling  with  flexure  overtones  and  weak  coupling 
with  face-shear  overtones  and,  in  addition,  weak  coupling  with  transverse 
thickness-shear  overtones.  As  for  the  interaction  of  the  3r<*  harmonic  thick¬ 
ness-shear  overtone  with  the  anharmonic  overtones  of  the  fundamental  thickness 
shear,  the  coupling  is  moderately  strong  at  small  a/b  (thick  plates  and  low- 
order  anharmonic  overtones)  and  diminishes  as  a/b  Increases  (thin  plates  and 
increasing  order  of  anharmonic  overtones). 

Finally,  the  minimum  absolute  values  of  the  slopes  of  the  segments  1^1 
are  much  larger  than  those  of  3^1  .  For  large  a/b  ,  the  ratio  of  those 
slopes  is  approximated  by  the  ratio  of  the  curvatures  of  branches  3^  and  lj 
at  za0  in  Fig.  1.  The  exact  values  of  those  curvatures,  in  the  three- 
dimensional  theory,  were  given  by  Eksteln  [7,  Eq.  56]: 
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<  *  [d2n/dz2]  *  k  +Ccot(rm/2  c^)  +  D  cot(mr/2  c|)  , 

n  z=0  c  i 

where 

k  =  (cu  +  A  +  B)  »  n  *  1  >3,5..., 

A  *  [(l+c12)cos  6  +  (c14  +  cg6)sin  0]2/(l-c2)  , 

B  *  C(cM  +  Cggjcos0  -  (l+c12)sin  6]2/(l-c3)  , 

C  *  4[(c2  +  c12)cos0  +  (c2c56  +  cu)]2/n2ir  c|(l-c2)2  , 

D  *  4[(c3  +  c12)sin  0  -  (c3cgg  +  c14)]2/n2iT  c|(l-c3)2  , 

c2»c3  =  ^2?  +  ^44  —  C(^22  -  c44)  +  4c24]^*}/2  , 
tan  0  *  c24/(c2  -  c44)  . 


For  the  present  case,  the  curvature  ratio  kj/k3  is  4.7  and  that  Is  the 
ratio  of  the  slopes. 

The  large  ratio  of  slopes  and  the  absence,  at  large  a/b  ,  of  strong 
coupling  with  all  overtones  except  those  of  flexure  (which,  at  such  high  over¬ 
tones,  have  very  small  amplitudes)  are  important  contributors  to  the  high 
stability  of  third  harmonic  overtone  resonators. 
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